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ABSTRACT 

The reception of FM carriers is characterized by a capture effect whereby the 
message of the dominant carrier is recovered when two or more FM carriers are 
present. This research uses a computer to form a type of average of the instan- 
taneous frequency of the receiver input. The results establish that averaging 
(smoothing or lowpass filtering) of the instantaneous frequency reveals the capture 
effect. The effect on capture of bandwidth and amplitude ratio for the case of two 
carriers is revealed. The results show that capture can occur when two carriers are 
separated by as little as 0.17 dB in amplitude. 
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I. INTRODUCTION 

Radio communication usage continues to expand because the military and 
commercial markets seek (1) personalized service such as paging and cellular 
applications, (2) ever greater exchanges of text and data (facsimile, computer 
networking), and (3) transfer of high quality signals (television). The messages are 
"traffic". The carriers are sine waves. 

There are two primary ways of carrying messages using sinusoids. The message 
can cause to vary (modulate) the amplitude or the argument (angle) of the sinusoid. 
These are called amplitude modulation (AM) and angle modulation (-^^M) 
respectively. 

Amplitude modulation can assume various forms such as double sideband with 
carrier, double sideband suppressed carrier and single sideband. Angle modulation 
can be either frequency modulation (FM) or phase modulation (PM). With analog 
messages, there is little difference between FM and PM. In fact, FM hardware 
would be used should PM be chosen. Consequently, PM with analog messages is, in 
practice, reduced to FM. [ Ref. 1: pp. 198—199 ] 

Frequency modulation exhibits several interesting characteristics not shared by 
AM. Two such characteristics are threshold effect and capture effect. The threshold 
effect relates to the quality (signal— to— noise ratio) of the demodulator output as a 
function of the quality of the received signal (demodulator input). Below 
"threshold" quality of the input, the output quality deteriorates rapidly. The 
capture effect relates to the ability of the demodulator to recover the message of the 
dominant carrier when two or more FM carriers are present. 
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This research is concerned with the capture effect of FM. Computer simulation 
is used to verify an idea concerning the reason for (origin of) the capture effect of 
FM demodulators. 

Background information on FM and the capture effect is presented in chapter II 
of this thesis. Chapter III describes the computer simulation procedure and 
program. Results of the simulation are contained in Chapter IV. Conclusions and 
recommendations are listed in chapter V. The computer program is listed in the 
Appendix. A list of references is provided. 
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II. BACKGROUND 



A. EQUATION FOR FM CARRIER 

An angle modulated sinusoid is expressed as 



s(0 = ficos[(^.(i)] (2.1) 

where the constant B is the carrier peak amplitude. A complete oscillation cycle 
occurs whenever the argument (pH) changes by 2:r radians. The instantaneous 
frequency of the angle modulated wave s{t) is defined by 



A 1 



( 2 . 2 ) 



Frequency modulation is that form of angle modulation for which the 
instantaneous frequency f^t) varies in accordance with the baseband signal or 
message voltage m{t). Typically, a linear relationship is chosen so that 



flt) = f^+ kfm{i) (2.3) 

where is the frequency of the unmodulated carrier and the constant kj represents 
the frequency sensitivity (in Hertz per volt) of the modulator. From Eqs. (2.2) and 
(2.3), the angle <j>{t) is 

(p{t) = 2jf t+ 27rkf f^m{T)dr (2.4) 

l C J J Q 



where, for convenience, the angle of the unmodulated carrier wave is assumed to be 
0 at i = 0. The frequency modulated wave is, therefore. 



s(^) = B cos[2x/^i + 2TikrJ^m{T)dT] 

J n 



(2.5) 
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B. CAPTURE EFFECT 

Consider the case of two FM carriers, si(<) and S 2 {t), present at the demodulator 
input and in the same frequency band. It is well known in practice that the output 
of the demodulator is the message of the dominant carrier. That is, if 

si(^) = Bi cos[27r / 1 H- 2’Kkr f^mi{T)dT] 

C 

S 2 {t) = B 2 cos[2;r(/ + €)t + 2nkf f^m 2 {r)dT] (2.6) 

c ]Jq 

then the output of the demodulator is proportional to 

when B\>B 2 
m 2 {t), when B\<B 2 

The carrier frequency offset c is any value such that both si(<) and S 2 (^) are in 
the same operating frequency band of the receiver /demodulator. This inherent 
ability of an FM demodulator to suppress weaker FM carriers is called capture 
effect. 

The FM demodulator can be a phase locked loop, a slope detector which 
includes Foster— Seeley and ratio detectors, or a pulse counting discriminator. 

C. ANALYSIS 

As shown by Eq. (2.3), the instantaneous frequency f^t) of a frequency 
modulated carrier is proportional to the message m{t) being carried. Consequently, 
it is claimed in the literature that the output voltage of a frequency demodulator is 
proportional to f^{t) — f^= [ Ref- 2: pp. 169—172 ]. This is indeed the case in 

practice when the demodulator input is a single frequency modulated carrier. 
However, it is not the case, because of the capture effect, when the demodulator 



4 



input is multiple frequency modulated carriers. 

The instantaneous frequency of the sum of two frequency modulated carriers 
can be derived as follows. Let 



v{t) = Bi COS(f>i{t) + B 2 COS02(O 
Using trigonometric identities, Eq. (2.7) can be written as 

v{t) = Ccos(f>z{t) 

where, 

+ B2'^ + 2BiB2COs[(f)i{t)-(l)2{t)] 

Sisin0i(/) + B2Sin(f)2{t) 

tan^3(<) = 

BiCOS<f>i{t) + B2COS(f)2{t) 

From Eq. (2.2), the instantaneous frequency f{t) is 









Bisin((>ifU+.52Sin((>2(^ 

i^lCO S<l>i{t)+ B 2 CO s4>2{t 



(2.7) 

( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 



_ (Bi0i'cos^i+.B2(^2'cos</i2 ) C+(.Bi<;i^i'sin<;;!»i+B2(?^2'sin(/>2)£> 

27r(C2+L>2) 

_ B\‘^<l>\+B2‘^(p2'+ [ BiB2COS{<l>i-<j>2)]{ <^l'~<^2') 

2 7t[5i2-|-52 2-+-251 52 COS ( ^l~02 )] 



( 2 . 11 ) 



where. 



C= Bi COS(f>i{t)+B2 COS(f>2{t) 
D= Bi sin(f>2{t)+ B 2 sin02(O 

/7 /7 

and <l>i, <f> 2 ' represent a.nd respectively. 
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As before, for an FM carrier 



= 2zfJ, + 2zkjJ^m.{T)dT 

= 27t/^ + 2nkjm.{t) i= 1,2 (2.12) 

where the carrier frequencies for both messages are assumed to be the same and 
m^{t) = i th message 

kj^ = frequency sensitivity of the frequency modulator in Hertz/volt 

Equations (2.8) through (2.12) can be reduced and interpreted for particular 
messages. For example, let the two messages be tones where 

m-{t) = D-cos2nf-{t) i= 1,2 (2.13) 

V I V 

For simplicity, let A = ^ = 1- Then, 



4>it) = 2-Kf t + 2jkf r^cos(27rf-r)dr 

Z C J J f\ z 



f . 

= 2nf t H sin27r/.f 

C t ’'Z 

h 



<t>' it) = 2-Kf + 2wk fCOs27rf-t i= 1,2 

Z C J z 

Using these values in Eq.(2.11) gives 



(2.14) 



where 



2 )+^ -X+Bi B2 Y- Z 

/.(() = ^ 

* A2+522+25i ^2^ 



(2.15) 



X = Bi2 cos27t Jit + B2^ cos2w fit 
Y = 2f^+kj-{cos2jft + cos2n fit) 

Z — cos[2nkj{ S\n2irf2t + 2 ^ Sin27r/U )] 
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Equation (2.15) is plotted in Figs. 2.1. a through 2.2.c along with Tri\{i) and 
M 2 {t) where /i = 1 Hz and /2 = 1-3 Hz. It can be seen that there are many 
impulselike frequency deviations from the message waveforms. The variable in Fig. 
2.1 is (3 ( modulation index for FM). The bandwidth of the modulated carrier is 
related to In general, bandwith increases linearly with (3 [ Ref. 3: pp. 194—197). 

There are two major factors that affect the graph of instantaneous frequency. 
First, the bandwidth of the modulated carrier affects the spacing in time between 
adjacent impulse— like frequency deviations. Figures 2.1. a through 2.1.C show' that as 
the bandwidth decreases, the spacing increases. Second, the difference in magnitude 
of the two modulated carriers affects the extent of the deviation. Figures 2.2. a 
through 2.2.C show that as the magnitude difference increases, the deviation extent 
decreases. 

However, these results do not show specific evidence of capture. Rather, they 
suggest that the instantaneous frequency needs to be averaged over a short time 
interval to reveal the capture effect. The same thing is true for the case of close (but 
different) carrier frequencies, as shown in Chapter HI. 

The idea that the capture effect is caused by the lowpass filtering operation of 
frequency demodulation led to this research. Since lowpass filtering can be thought 
of as averaging (or smoothing), it was decided to investigate the capture effect by 
averaging the instantaneous frequency using a computer. 

So, in hardware, the averaging is done by the lowpass filter (LPF). A LPF is 
part of every frequency demodulator (frequency— to— voltage or voltage— to-frequency 
converter). Consequently, the output of a frequency demodulator is a lowpass 
filtered version of instantaneous frequency deviation A f{<). 
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Again, lowpass filtering can be considered as a smoothing operation on the 
applied voltage. In this research, we used computer simulation to investigate the 
effect of smoothing of to reveal the presence of capture. Smoothing is created 
by averaging time intervals between adjacent zero crossings of the FM carrier. Since 
this interval (Ti) represents the time required for tt radians excursion of the 
argument of the cosine carrier, and since tt radians is taken as a half cycle, then two 
such intervals ( Ti + T2) become a full cycle and is the number of cycles per 

second or the instantaneous frequency for Ti + T2 seconds. Averaging N such 
groupings provides a smoothed value of or short-term average of the instantaneous 
frequency over that expanded time interval. 

A computer program was used to approximate Ti and T2 for every interval 
between zeros of v{t). These values are then averaged and the results are plotted. 
The program is explained in Chapter III. 
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0 0.2 0.4 0.6 0.8 1 
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Frequency 




Fig. 2.1 Instantaneous Frequency for Tone Messages 
where a^ja^ = 1.2, = 200 Hz, ^ = 1 Hz, = 1.3 Hz 
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Frequency 




T ime 

a. fli/(72 = 1 




Time 

b. Oi/c2=l.l 

) 
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Frequency 




Fig. 2.2 Instantaneous Frequency for Tone Messages 
where P = 10, f = 200 Hz, /i = 1 Hz, = 1.3 Hz 

C 
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m. SIMULATION PROCEDURE 

The FM demodulator, including LPF, is simulated by a computer program so 
that critical parameters, such as signal amplitudes, bandwidth of the FM system, 
message forms, and interval of averaging, can be easily changed. Since the sampled 
version of the demodulator input forms a sequence of numbers, it is convenient to 
consider it as a row (or column) vector. Once that is done, the rest is simple 
manipulation of the row (or column) vector. 

A. WAVEFORM GENERATION 

The desired waveform v{t) can be generated using a simple routine 
(WAVEGEN). In this routine, the carrier amplitudes, values of kj, carrier frequ- 
encies and message frequencies are input parameters which can be changed. This 
resulting wave is sampled at a rate which will provide the desired accuracy in 
measuring the time between adjacent zeros of v{t). In this work, a carrier frequency 
of 200 Hz and sampling rate of 12 kHz are used. 

The sampled values contain no information about frequency. Rather, the 
number of samples between zero crossings are related to frequency because the 
number of samples between zero crossings itself represents a half cycle oscillation of 
argument of v{t) and its reciprocal can be interpreted as twice the frequency of that 
period. Therefore, the routine notes the polarity of each sample so that the number 
of samples between zero crossings (polarity changes) can be easily counted. This 
result is then used as the input for the next routine. 
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B. COUNTING 

The number of samples between zero crossings (polarity changes) multiplied by 
the fixed sampling interval (A) is considered to be taken as half the instantaneous 
period of that interval. So, a sequence of integers which consists of the number of 
samples between zero crossings need only be generated. The i th entry for this 
sequence is denoted by , j = 1 , 2 , 3, - • • as shown in Fig. 3.1. The totality of 
entries forms a sequence of integers. 



Furthermore, T- itself is a time index. Since there are i— 1 entries before the 

i 



the interval T- is selected for the time index of interval T - . For example, the time 

i 



Since the first and last interval may contain a number of samples that does not 
correspond to a complete interval between adjacent zeros, these two entries are not 
used in the averaging routine. 




Fig. 3.1 Formation of sequence 




index of T 3 is 7\ + T 2 + ^^ 3 .) 
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C. AVERAGING 



Given the sequence of integers, any average Tj^ over N intervals is formed as 

k+N 
S T. 
k * 

T, = k = any positive integer (3.1) 

A 

T- = i th interval 

i 

N = number of intervals averaged 



Since frequency is taken as the reciprocal of the period, the averaged frequency, 
denoted by , is given by 




1 

2T^X 



(3,2) 



To obtain more realistic (smoother) results, a moving average is used rather 
than fixed time interval averaging. To illustrate this, assume a sequence of intervals 
as follows 



Ti T2 Ts T4 n 

Then, if the average is done over 3 intervals, the averaged version of the sequence 
will have the following entries 

rp _ T2+ Tz 

rp _ T2+TZ+T4 

I2- 3 




^k '^k+\ '^k+2 



(3.3) 
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Again from the relationship between period and frequency (Eq.(3.2)), averaged 

version of frequency will form the sequence 

l_ 1 ]_ _ 

2T X 2 T 2 X 2T^X 

In this case, the time indices for the averaged version of intervals are 

= time index for T\ 

Ti_|_ 72+-|3+Z4 = time index for T 2 

Ti+ T 2 + ^ 2 ~*~ = time index for and so on. 

Generally, time indices for the moving average over A intervals are found as 
i — 1 i-\- N-\ 

AA + 5 kh A (time index for ) (3.4) 

Since the values of each T\ are not the same, the averaged frequency values w’hen 
plotted will not be equally spaced. 
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IV. SIMULATION RESULTS 

With an average of 2 intervals {N = 2), the resulting plot is very similar to that 
of the analytical one. That is, when almost no averaging is done, the resulting plot 
should yield that of the analytical one. This is another way of verifying that the 
simulation program is running properly. As the number A of average intervals 
increases, the impulse— like frequency deviation becomes smaller and the overall 
frequency deviation approaches that of the stronger signal. 

Figures 4.1. a through 4.5.C are the simulation results. The results are presented 
as averaged instantaneous frequency versus time. In all cases, the message mi(t) 
corresponds to that of the dominant carrier. The following variables are used to see 
the effect of system parameters. 

aj/fl 2 = ratio of carrier amplitudes (1 through 1.5) 

0 = modulation index (1 for narrowband, 10 for wideband, 50 
for very wideband FM) 

^ = message 1 frequency (1 Hz for sinusoidal message) 
fi = message 2 frequency (1.3 Hz for sinusoidal message) 
fc = carrier frequency (180 Hz through 200 Hz) 

N = number of intervals averaged (2, 20, 40, 60, 80) 

Various messages were used in the simulation. Figures 4.1. a through 4.1.C show 
the results when the messages are linear intrapulse FM (up CHIRP and down 
CHIRP). Figures 4.2.a and 4.2.b are for exponentially decaying and growing 
sinusoidal messages. Sinusoidal messages are used extensively to show the results for 
various correlations of the listed variables. 
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Even though there may be other factors which influence the capture effect, only 
the following four factors are investigated in this research. 

A. EFFECT OF AVERAGING TIME) 

Figures 4. 3. a through 4.3. d shows the effect of N. The other factors ( 01 / 02 , (i, Ji, 
fc) remain fixed. This factor is essentially the cause of the capture effect. As N 
increases, the instantaneous frequency increasingly resembles the message of the 
stronger signal. 

B. EFFECT OF SIGNAL AMPLITUDE RATIO 

Figures 4.4.a through 4.4. d show the effect of signal amplitude difference. While 
the other factors remain fixed, 01/ 02 varies from 1.01 to 2.0. The effect of signal 
amplitude difference is very interesting. Even with very small signal amplitude 
difference, the system still locks onto the stronger signal. For numerical comparison, 
the 0.8 dB difference case (Fig. 4.4.b) and 1.6 dB case (Fig. 4.4. c) show little 
difference. In the case of very small amplitude difference, more averaging is needed 
to see the capture effect. Figure 4.4. a indicates capture for a carrier separation of 
0.17 dB (ai/fl2 = 1.02). 

C. EFFECT OF MODULATION INDEX 

Modulation indicies are changed from 1 to 50 in Figs. 4.5. a through 4.5.c. 
Capture occurs for both narrowband FM and wideband FM. Since the narrower 
band FM shows fewer impulse— like frequency deviations, each of longer duration 
than those of wider band FM, narrow band FM needs to be averaged longer (see 
Fig. 2.1). 
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D. EFFECT OF CARRIER FREQUENCY DIFFERENCE 

Figures 4. 6. a and 4.6.b show the effect of carrier frequency difference on 
capture. It appears that capture becomes more complete as the carrier frequencies 
separate. Obviously, the two carriers must lie in the passband of the receiver 
(intermediate frequency amplifier) if they are to be present at the input to the 
demodulator. 
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Frequency 




a. /V= 2 




b. N=20 
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Frequency 




c. /V=40 

Fig. 4.1 Average insLantaneoiis Frequency Deviation for CHIRPS 

where 01/02 = 1.5 
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b. /V=20 

Fig. 4.2 Average Inslanlaneous Frequency Deviation 
for Decaying and Growing Sinusoidal Messages where ai/(i 2 = 1.2 
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Frequency 




T ime 

a. yv= 2 




b. yv=20 
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Frequency Frequency 



220 





d. A^=80 



Fig. 4.3 Average Instantaneous Frequency Deviation for Tone Messages 
where asja^ = 1.2, /? = 10> 4 = 200 Hz, fi = 1 Hz, /i = 1.3 Hz 
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Frequency 





( 1 . ffi / 02 = 2.0 

Fig. 4.4 Average Inslanlaneous Frequency Deviation for Tone Messages 
where /? = 10, /^ = 200 llz, /, = 1 Hz, ^ = 1.3 Hz, A = 60 
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a. P = I 




T ime 



b. /?= 10 
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Frequency 



260 




c. /?=50 

Fig. 4.5 Average Instantaneous Frequency Deviation for Tone Messages 
where 01/02 = 1.2, = 200 IIz, /i = 1 Hz, ^ = 1.3 Hz, = 60 
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Frequency 





b. /cl = 190117., /c 2 = 200Hz 



Fig. 4.6 Average Inslanlaneous Frequency Deviation for Tone Messages 
where m/oj = 1-2, /, = 1 IIz, /2=1.3 IIz, /? = 10, /V= 60 
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V. CONCLUSIONS AND RECOMMENDATIONS 



The simulation results establish that the lowpass filtering portion of frequency 
demodulation accounts for the capture effect of FM receivers. Without filtering, the 
instantaneous frequency deviation of the sum of two carriers has an impulse— like 
history which does not consistently favor the message of the dominant carrier. With 
averaging or smoothing (lowpass filtering) of the instantaneous frequency deviation, 
capture occurs. 

The simulation results indicate capture occurs for both narrowband and 
wideband FM. Capture is evident when the carriers are separated in amplitude by 
less than 0.2 dB. The effect of averaging time is revealed by the simulation output. 
In general, longer averaging times produce smoother approximation to the message 
of the dominant carrier, as expected, and also show the capture effect more 
convincingly. 

It is recommended that the capture effect be verified using an operating 
experimental system in which system parameters can be controlled and accurately 
measured. 
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APPENDIX 



Program List 

This appendix contains the program used to calculate f- defined in the thesis. 

V 

The program actually determines T^ and as given by Eq. (3.1). 

The input equation is Si(t) + S 2 (t) as defined by Eq.2.6. This equation is 
sampled 12,000 times in the one second simulation interval. The interval between 
samples is, then, 1/12,000 = 83.3 fisec = A. The polarity of each sample value is 
noted. The number of consecutive sample values of like polarity in interval i is . 

The program stores each value of T. and forms T^, from Eq. (3.1) for various 
assigned values of N. Then ^ is calculated according to Eq. (3.2). Finally, ^is 
plotted for each time index of Eq. (3.4). 

The program contents are following: 

MAIN : main routine for the simulation of capture effect 

CHIRP GEN : subfunction routine for generation of CHIRPs 

WAVEGEN : subfunction routine for generation of sinusoidal waves 

COUNT : subfunction routine for counting and averaging 

SAMPLE : subfunction routine for display 

The program was written in APL and run on IBM PC/AT with 3 Mb of memory. 
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A. MAIN 



V Y 4-MAIN A;A; Y ;TMPSMP;TSLICE 

[1] n 

[2] n This is main routine for the capture effect 

[3] n Simulation. For the clality of logical flow, 

[4] n each simulation step is divided into their 

[5] n own routines. 

[6] n The output will be, 

[7] n Ts time index for fl and f2 

t8] n TIME: time index for aceraged frequency 
[93 n (not equally spaced) 

[103 n F1,F2: individual instantaneous 

[113 n frequency 

[123 n RES: averaged version of frequency 

[133 n 

[143 

[153 Q4- ' Genera ting waves....' 

[ 163 

[173 n Generating arbitrally sinusoidal waves 
[183 

[19 3 TMPSMP4-WAVEGEN A 
[203 ' Counting . . . . ' 

[213 

[223 n Count the number samples between 

[233 n zero crossings and form a resulting sequency 

[24 3 

[25 3 Y4-C0UNT TMPSMP 
[263 D4-' Sampling. ... ' 

[273 ^ ^ 

[283 n Since the numbers of elements of sampled time and 
[293 n message are too many for plotting purpose, it is 
[303 n sampled without missing any important information. 
[313 

[32 3 T4-S AMPLE T 
[ 33 3 F14-SAMPLE Fl 
[ 34 3 F24-S AMPLE F2 

[ 35 3 U<-'Done.' 

V 
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B. CHIRPGEN 



V V4-CHIRPGEN A1;DI0;A1;A2;PH1;PH2;V;TMP1;TMP 

[ 1 ] 

[2] fi This routine generates two CHIRP signals, namely up 

[3] 0 and down CHIRP to test the program. 

[4] R Equations for this routine are, 

[5] n F1=100T+150 for up CHIRP 

[6] n F2=250-100T for down CHIRP 

[7] fi where OsT^l 
C8] 

[9] DI04-0 

CIO] A24-1 

[11] TSLICE^-r (1. 3*0. 5)xl0000 

[12] T^(lTSLICE)-^TSLICE 

[13] D«-50 

[14] TMP^(02)xD 

[15] TMP4-TMPXT 

[16] TMP4-TMPXT 

[17] TMP1^(O2)xC200-D) 

[18] TMP14-TMP1XT 

[19] PH1*-TMP+TMP1 

[20] TMP<-CO2)x(200 + D) 

[21] TMP4-TMPXT 

[22] TMPl<-C02)xD 

[23] TMPl^-TMPlxT 

[24] TMP14-TMP1XT 

[25] PH2+-TMP-TMP1 

[26] TMP4-20PH1 

[27] TMP14-20PH2 ^ 

[28] TMP<-A1xTMP 

[29] TMP1^A2xTMP1 

[30] V4-TMP+TMP1 

[31] V<-xV 

[32] F1^2xDxT 

[33] FK-F1+200-D 

[34] F2+-200 + D 

[35] TMPf-2xD 

[36] TMP^TMPxT 

[37] F2«-F2-TMP 
V 
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C. WAVEGEN 

V V«-WAVEGEN A1 ; EIO ; A1 ; A2 ; V ; TMPl ; TMP ; VI ; V2 ; KFl ; KF2 
t 1 ] ; FCl ; FC2 ; FMl ; FM2 ; PHI ; PH2 ; PI ; P2 

t2] 0 

t3] 0 Routine to generate samples of sinusoid waveforms 

[4] fi 

t5] DIO<-0 

t6] 

t7] n Second wave amplitude. This routine accepts the first 

[8] n signal amplitude as an input argument 

[9] 

[10] A2<-1 

[11] P2^-o2 

[12] TSLICE<-r (1. 5*0. 5)xl0000 

[13] T*-C VTSLICE)-^TSLICE 

[14] 

[15] 0 Carrier frequency, frequency sensitivity and message 

[16] fi frequency 

[17] 

[18] FC14-200 O FC24-200 

[19] KFK-IO O KF2<-10 

[20] FM1*-1 O FM24-1.3 

[21] PHl«-0 O PH2<-0 

[22] TMP*-P2xFMl 

[23] TMP4-TMPXT 

[24] TMP4-TMP+PH1 

[25] TMP<-10TMP 

[26] TMP<-TMP-( lOPHl) 

[27] TMP«-TMPxKFl^FMl 

[28] TMP14-P2XFC1XT \ 

[29] TMP+-TMP+TMP1 

[30] VK-20TMP 

[31] V14-V1XA1 

[32] TMP+-P2XFM2 

[33] TMP4-TMPXT 

[34] TMP*-TMP+PH2 

[35] TMP*-lOTMP 

[36] TMP«-TMP-(10PH2) 

[37] TMP*-TMPxKF2^FM2 

[38] TMPl*-P2xFC2 

[39] TMPK-TMPlxT 

[40] TMP<-TMP+TMP1 

[41] V2<-20TMP 

[42] V2*-A2xV2 

[43] V<-x(Vl+V2) 

[44] TMP4-P2XFM1 

[45] TMP*-TMPxT 

[46] TMP4-TMP+PH1 

[47] TMP*-KF1XC0S(TMP) 

[48] FK-FCl+TMP 

[49] TMP*-P2xFM2 

[50] TMP*-TMPxT 

[51] TMP4-TMP+PH2 

[52] TMP*-KF2xCOS(TMP) 

[53] F2«-FC2 + TMP 
V 
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D. COUNT 

V Y4-COUNT V;C;COUNTER;INIT;I;LIMIT;D;Y;TMP 

[1] ;EXTIME1;EXTIME2 ;INTV 

[ 2 ] 

[3] nRoutine to form half period sequence and average it 

[4] 

[5] fi 

C6] R Set the initial matricies and values 

C7] n 

t8] 

C9] COUNTER<-0 
CIO] INIT<-VCDIO] 

C113 I4-DI0 

C12] LIMIT<-(pV)+DIO-l 
[13] n 

C14] fi Routine to count number of samples between 

[15] R Zero-crossings 

[16] R 

[17] LOOPl: 

[18] ->(V[I]=INIT)/L00P2 

[19] C<-C, COUNTER 

[20] INIT<-V[I] 

[21] C0UNTER<-1 

[22] ->LOOP3 . 

[23] LOOP2: 

[24] COUNTER<-COUNTER+l 

[25] LOOP3: 

[26] I<-I + l 

[27] ->(ISLIMIT)/L00P1 

[28] C<-C, COUNTER i 

[29] R 

[30] R Drop the incomplete cycle at the beginning and the end 

[31] R 

[32] EXTIME1«-C[DIO]x(14^TSLICE) 

[33] C<-liC O C4-”11C 

[34] D<-0 

[35] TIME<-^ 

[36] R 

[37] R Average with the value of INTV 

[38] R 

[39] INTV4-60 

[40] I<-DIO 

[41] LIMIT<-CpC)-(UIO + l + INTV) 

[42] R 

[43] R Calculate extra time due to averaging 

[44] R 

[45] EXTIME2<-( + /INTV1C)^TSLICE 

[46] R 

[47] R Calculate the time index which is not equally placed 

[48] R 

[49] LOOP4: 

[50] TMP<-+/INTVTC 

[51] D<-D,TMP^INTV 

[52] TIME<-TIME,C[DIO] 

[53] C<-liC 

[54] I<-I + l 
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t55] ->(I<LIMIT)/L00P4 

[56] TIME<-( ( + \TIME)■^TSLICE)+EXTIMEl+(EXTIME2-^2) 

[57] n 

[58] n Calculate averaged frequency for the time indicies 

[59] fi 

[60] Y<-(1^2xD)xTSLICE 
V 
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E. SAMPLE 



V Y<-SAMPLE X;DIO;Y;SIZ;LIMIT;I;X 

[ 1 ] 

[2] R Routine to sample for plotting (number of samples = 60 

[3] 

[4] D10<-0 

[5] Y<-^ 

[6] SIZ<-L ( (pX)^600) 

[7] LIMIT+-(pX)+DIO-l 
[83 I<-DI0 

[93 LABELl: 

[103 Y<-Y,X[lxSIZ3 
[113 I<-I + l 

[123 -♦( (IXSIZ)SLIMIT) /LABELl 

V 
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